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I.  INTRODUCTION 


Application  requiring  combustion,  power  generation  or  heat 
transfer  on  land,  on  the  air  and  in  space  by  methods  over  than 
direct  conversion  (photovoltaic)  will  require  new  concepts 
materials  to  obtain  optimum  performance  under  desired  high 
operating  temperatures.  Traditional  materials  such  as  zircaloy, 
stainless  steel,  precipitation  and  dispersion  strengthened  steels 
and  on  some  exotic  applications,  Nb  and  V  steels,  do  not  have  the 
necessary  combined  ultrahigh  temperature  strength,  irradiation 
resistance  and  strength  to  weight  ratios  that  are  desired. 
Operating  temperatures  in  excess  of  1500  degrees  C,  where  most  of 
the  above  traditional  materials  melt,  are  desired. 

Carbon  based  materials  (C  materials)  are  a  superb  candidate 
for  development.  C  fiber-C  matrix  materials  have  been  highly 
successful  on  the  space  shuttle  where  they  have  been  used  for 
components  ranging  from  cargo  bay  doors  to  high  temperature  motor 
casings.  Such  materials,  however,  are  not  suited  for  sustained 
high  temperature  application  or  to  cases  which  require  exposure 
to  nuclear  radiation  or  long  term  exposure  to  fluid  or  gas 
pressure.  Testing  in  high  temperature  gas  cooled  reactors  has 
identified  two  problems  with  C  materials. 


The  first  is  the  permeability  of  the  C  material  to 
radioactive  fission  gases  [1,2].  Despite  the  development  of 
ultrahigh  density  C  materials,  special  C  coatings,  and 
differential  pressurization  of  components,  leakage  of  C  clad  fuel 
occurs  after  very  short  neutron  exposures (^t  <  1  x  1021  n/cm2 ) . 
Heat  pipe  components,  fabricated  from  pure  C  or  C  composite 
materials,  would  therefore  be  unable  to  maintain  proper  partial 
pressure  for  operation.  The  second  problem  results  from 
volumetric  change  (desification)  that  occur  in  C  materials  upon 
neutron  irradiation  [1] .  These  changes  are  highly  orientation 
dependent.  The  C  fibers,  which  contain  material  with  the  basal 
planes  aligned  perpendicular  to  the  stress  axes,  undergo  a 
substantially  different  volumetric  change  than  the  matrix 
material,  leading  to  the  generation  of  high  internal  stress, 
cracking  and  ultimate  failure. 

A  class  of  materials  which  promise  to  overcome  the  above 
described  problems  is  ordered  C-metal  alloys.  Ordered  alloys 
offer  many  potential  advantages  over  conventional  alloys  at 
elevated  temperatures  [2-5].  The  atomic  ordering  produces  a 
pronounced  increase  in  work  hardening,  [5-9]  improves  the  fatigue 
resistance,  [10]  and  retards,  because  of  stronger  binding  and 
closer  packing  of  atoms,  thermally  activated  processes  such  as 
creep  and  grain  growth  [11].  In  addition,  the  strength  of 
ordered  alloys  is  less  sensitive  to  temperature  than  conventional 
disordered  alloys.  In  fact  some  alloys  show  an  increase  rather 
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Chan  a  decrease  in  strength  with  increasing  temperature  up  to  the 
critical  ordering  point  [12].  The  critical  ordering  point  is  the 
temperature  at  which  the  material  reverts  from  a  defined  periodic 
arrangement  of  alloying  elements  to  the  random  arrangment  found 
in  conventional  materials.  Electron,  ion  bombardment  and  recent 
neutron  irradiation  results  show  the  long  range  ordered  (LRO) 
alloys  to  be  highly  swelling  resistant  [13-15]. 

Despite  the  above  advantages ,  the  LRO  alloys  have  seen  only 
limited  application  because  of  a  lack  of  ductility  associated 
with  the  ordered  state  [16-18].  Recently  Liu  and  his  coworkers 
at  ORNL  have  suceeded  in  producing  (Fe,Ni)3V  long  range  ordered 
alloys  which  show  tensile  elongations  greater  than  30%  [19-21]. 
This  has  been  achieved  by  controlling  the  ordered  lattice 
structure  though  use  of  the  e/a  ratio  or  average  electron  density 
per  atom  outside  the  inert  gas  shell.  In  other  words  by  creating 
off  stoichiometric  compounds . 

Ordered  C -metal  alloys  represent  a  new  approach  which 
combines  the  desirable  qualities  of  C  (high  temperature 
resistance,  low  weight,  low  vappor  pressure)  with  the  radiation 
and  permeability  resistance  of  ordered  materials.  These  alloys, 
containing  38-60  atomic  percent  carbon,  all  exhibit  melting 
temperature  in  excess  of  2500  degrees  C  [22,23].  By  careful 
control  of  the  alloy  composition,  desirable  ordered  structures 
can  be  maintained  at  temperatures  in  excess  of  0.7  Tm  (approx. 
1800  degrees  C) .  The  alloys  will  be  relatively  lightweight,  yet 


resistant  to  irradiation  damage.  Control  of  the  ordered 
structure  by  modification  of  the  e/a  ratio  will  produce 
compositions  with  good  ductility,  strength,  creep  and  fatigue 
properties.  In  addition,  the  C-metal  alloys  will  eliminate  the 
permeability  problem  associated  with  pure  C  materials. 

The  C-Ti  system  based  upon  its  ability  to  obtained  an 
ordered  structure  over  a  range  of  compositions  and  its  ultrahigh 
melting  temperature  (3080  degrees  C) .  The  objective  of  this 
study  is  to  investigate  the  structure-property  relationships  of 
0-40  weight  percent(w/o)  carbon  of  the  C-Ti  system.  The 
experiments  concentrate  on  an  investigation  of  the  effect  of  C 
content  on  the  mechanical  properties.  These  studies  provide  a 
theoretical  and  experimental  basis  for  understanding  the 
properties  of  C  base  alloys. 
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II.  THEORY  OF  ORDERED  ALLOYS 


There  are  a  number  of  theoretical  models  available  which 
allow  the  calculation  of  the  critical  ordering  temperature, 
maximum  degree  of  order  and  specific  heat  as  a  function  of 
temperature  for  ordered  alloy  systems.  It  is  the  objective  of 
this  research  to  apply  these  theories  to  the  C-Ti  alloy  system. 
The  results  of  the  above  computations  will  provide  a  theoretical 
basis  upon  which  the  most  promising  alloy  systems  and  ranges  of 
compositions  can  be  chosen  for  melting  and  experimentation.  It 
is  hoped  that  before  alloying  additions  are  made  that  we  will  be 
able  to  extend  these  theoretical  calculations  to  ternary  systems 
[24]. 


I*. 


Unlike  conventional  alloys ,  Which  from  a  random  mixture  of 
atoms  on  lattice  sites,  long  range  ordered  (LRO)  alloys  exhibit 
periodic  arrays  which  form  an  ordered  crystal  structure.  In 
general,  unlike  atoms  tend  to  be  closer  together  and  like  atoms 
further  apart.  The  tendency  to  order  increases  with  decreasing 
temperature,  therefore,  the  ordered  arrangement  of  atoms  must  be 
more  stable  or  energetically  favorable  at  lower  temperature  than 
a  random  distribution.  For  this  to  be  ture,  the  attraction 
between  unlike  atoms  stronger  than  the  attraction  between  like 
atoms.  As  the  temperature  of  the  alloy  increases,  the  energy  of 
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thermal  vibration  increases  and  when  this  exceeds  the  so  called 


ordering  energy,  the  atoms  have  sufficient  energy  to  redistribute 
themselves  into  a  disordered  or  random  mixture.  It  is  at  this 
point  that  the  disordered  structure  becomes  the  stable 
configuration.  Thus,  we  expect  a  transition  temperature  at  which 
the  ordered  phase  becomes  unstable.  This  transition  temperature, 
the  critical  temperature,  Tc,  is  analogous  to  the  Curie 
temperature  in  ferromagnetism. 


Critical  Ordering  Temperature 
To  predict  the  critical  ordering  temperature  for  C-Ti 
alloys,  a  review  of  various  theoretical  models  of  ordering  was 
undertaken.  This  review  is  summarized  in  what  follows.  Above 
the  critical  temperature  the  alloy  is  completely  disordered. 

Once  the  critical  temperature  is  reached,  atomic  ordering  begins 
and  at  each  temperature  below  the  critical  temperature  (after 
suitable  time)  an  equilibrium  degree  of  order  will  be  reached. 
Thus,  at  the  stoichiometric  composition  of  the  alloy,  the 
equilibrium  degree  of  order  at  high  temperatures  will  be  zero  and 
at  0  degree  kelvin  the  alloy  will  be  completely  ordered,  alloys 
which  are  not  stoichiometric  will  not  reach  perfect  order,  but 
will  have  a  characteristic  maximum  degree  of  order  depending  on 
the  composition. 

the  first  suggestion  that  some  alloys  can  exhibit  an  ordered 
structure,  or  supper-lattice  was  put  forth  in  1919  by  Tammann 


[25]  after  his  observations  in  the  Cu-Au  system.  Since  that  time 
a  number  of  efforts  have  been  made  to  theoretically  describe  and 
predict  the  cooperative  phenomena  of  ordering  [26]. 

The  first  successful  theory  of  the  stability  of 
superlattices,  as  a  function  of  temperature,  was  developed  by 
Bragg  and  Williams  [27].  Later,  Bethe  [28]  developed  a  more 
refined  method,  which  was  extended  by  several  authors.  Neither 
one  of  these  methods  were  very  clear  about  the  assumptions  which 
were  made;  Kirkwood  [29]  remedied  this  by  developing  a  theory 
based  entirely  on  the  standard  concepts  of  statistical 
mechanics.  The  quasichemical  method  was  developed  also  at  this 
time  by  Guggenheim  [30]  and  was  later  shown  to  be  equivalent  to 
Bethe ' s  method  [ 31 ] . 

The  more  theories,  including:  the  Cluster-variation  method 
developed  by  Takagi  [32],  Kikuchi  [33],  Yang,  Li  [34]  and  Hill 
[35],  the  method  of  concentration  waves  [11],  and  a  band  theory 
model  [37].  In  addition  to  these,  numerical  methods  based  on  the 
random  walk  or  Monte  Carlo  Method  [38]  have  been  developed. 

Except  for  band  theory,  all  of  the  above  models  require  the 
interaction  or  ordering  energies  of  the  alloy, 

Vl  -  V*A  +  V*B  -  V*B  . (1) 

where  the  superscript  i-1  denotes  nearest  neighbors,  i-2  denotes 
next  nearest  neighbors,  etc.,  and  the  subscripts  denote  which 
pair  of  atoms  the  energy  Vap  represents.  Unfortunately,  this 
energy  parameter  is  hard  to  come  by,  This  will  be  discussed  later 
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Bragg-Williams ,  Kirkwood,  Quasichemical  Methods 
To  specify  the  degree  of  order  of  the  arrangement,  Bragg  and 
Williams  introduced  the  long  range  order  parameter  s  in  their 
original  paper  [27].  When  perfect  order  is  achieved  in  a  binary 
alloy  of  components  A  and  B,  the  lattice  sites  occupied  by  A 
atoms  are  called  a  sites  and  those  occupied  by  B  atoms  are  called 
p  sites.  Let  N  be  the  total  Number  of  atoms  (A  and  B)  and  thus 
the  total  number  of  lattice  sites,  (i.e.  no  vacancies);  and  FA 
and  Fb  (where  FB  -  1-  FA )  be  the  fraction  of  A  and  B  atoms  in  the 
alloy  respectively.  Also  let  the  fraction  of  a  sites  occupied  by 
A  atoms  be  denoted  by  ra  (right  atoms),  the  fraction  of  a  sites 
occupied  by  B  atoms  be  wa,  where  wa  -  l-ra  and  similarly  define 
rp  and  wp  -  1-r^  represent  the  rightly  and  wrongly  occupied  p 
sites.  The  long  range  order  parameter  s  may  be  defined  as 


Thus,  s  is  defined  as  unity  for  perfect  order  and  vanishes  in  the 
completely  disordered  state. 

Now,  let  Q  be  the  total  number  of  all  types  of  nearest 
neighbor  pairs  in  the  lattice  and  each  atoms  be  surrounded  by  z 
nearest  neighbors,  thus 

Q  -  (Z/2)  N  .  ---(3) 

Also  let  the  number  of  pairs,  which  are  AA,  BB  and  AB  be  denoted 
by  Qaa ,  QBB ,  QAB  respectively.  Consider  an  alloy  with  FAN  a 
sites  and  FgN  p  sites.  Then  according  to  equation  (2)  the 
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average  number  of  A  and  B  atom  on  a  and  0  sices  are  as  follows 
[26]: 

A  atoms  on  a  sices:  raFAN  -  (FA  +  FBs)FAN, 

A  atoms  on  0  sites:  w^FBN  -  (1  -  s)FAFBN  , 

B  atoms  on  a  sites:  wQFAN  -  (1  -  s)FBFAN  , 

B  atoms  on  0  sites:  r^FBN  -  (FB  +  FAs)FBN, 

The  number  g(s)  of  distinguishable  arrangements  or  configurations 
for  a  given  state  s  is  then 

[FaS]  [FbN] 

g(a) - - 

[ Fa  +Fb s ) FA N ]  !  [  (1-s)FaFbN]  !  [  (l,s)FAF„N]  !  [  (FB+FAs)FBN]  ! 

. (5) 

Assuming,  that  the  lattice  vibrations  are  independent  of  the 
configuration  of  the  atoms  in  the  lattice,  the  partition  function 
in  the  canonical  ensemble  may  be  written  as 

z(T, s)  -  S  e*wr/KT  . (6) 

where  T  is  the  absolute  temperature,  s  is  the  long  range  order,  K 
is  the  Boltzmann  constant  and  Wr  is  the  configurational  energy  of 
the  state  r .  The  summation  extends  over  all  the  state  r  of  a 
given  order  s. 

Our  knowledge  of  Wr  is  scant  at  best  at  this  time  and  is 
approximated  here  as  consisting  mainly  of  the  interaction 
energies  of  pairs  of  nearest  neighbors.  These  energies  as 
identified  in  eq.(l)  are  VAA ,  VBB ,  VAB ,  corresponding  to  nearest 
neighbor  pairs  QAA ,  QgB ,  QAB  respectively.  Therefore  the 
configurational  energy  is  roughly  given  by 


for  any  specified  configuration  r.  The  configurational  free 
energy  F(s)  is  given  by 


(8) 


F(s)  -  -KTln(Z(T , s) )  , 

and  once  this  is  minimized  the  equilibrium  value  of  s  will  be 

determined.  The  corresponding  energy  of  the  alloy  is 

ZWre-“r0 
r  r 

E(s)  -  W(s)  -  .  . (9). 

2  e-“r0 

T 

where  p  is  1/KT. 

Now,  using  the  conventional  procedures  of  statistical 
mechanics,  we  define  U(s)  by  relation 

2  e*uT0  -  g(s)e*u (• . (10) 

Here  g(s)  is  the  same  as  in  eq.  (5).  Now,  eq.  (8)  can  be  written 
as 

F(s)  -  U(s)  -  KTln(g(s))  ,  . (11) 

and  the  energy  of  the  alloy  can  be  written  as 


8  U(s) 

E(s) . . (12) 

3(1/T)  T 


The  above  equations  for  the  free  energy  and  internal  energy 
define  the  thermodynamics  of  the  system  [26],  All  that  is  left 
to  do  is  to  calculate  the  partition  function.  This  however  has 
proved  illusive  and  various  approximations  have  been  made,  the 
first  of  which  was  the  Bragg-Williams  approximation  [27],  Their 
approximation  consists  of  replacing  the  true  value  of  Wr  for  each 


configuration  of  s  by  the  simple  average  <Wr>  of  all  staets 
belonging  to  a  given  s, 

<wr>  -  -  {<Q AA>VAA  +  <088^88  +  <QaB>VAbJ  >  --(13) 

where  <Q*a>,  <Qes>  anc*  <Qab>  are  t^'e  sample  averages  of  the 
number  of  pairs  of  each  kind  in  each  state  r  of  s.  Consider 
<Qaa>,  AA  pairs  may  result  from  two  arrangements:  (1)  one  A  atom 
on  an  a  site  and  other  one  on  a  0  site,  and  (2)  both  atoms  on  0 
sites.  On  the  average  there  are  FAN(FA  +  Fas)a  sites  occupied  by 
A  atoms.  There  are  then  zNF|(1-s)(Fa  +  FBs)  AA  pairs  of  type 
(1) .  The  number  of  AA  pairs  with  both  atoms  on  0  sites  is  1/2 
zNFj (Fb-Fa)(1-s)2  [26].  Therefore, 

<Qaa>  -  zNF|(1-s)(Fa+Fbs)  +  1/2  zNFf  (F„-FA)(l-s)2 

-  1/2  zN(F|-F|s2)  . (14) 


and  similarly, 

<Qb8>  -  1/2  zN(F§  -  Fts2) 

<Qab>  -  1/2  zN(2FaFb  +  2Fa s2  )  . 
Substituting  equations  (14)  into  (13)  leads  to 

<Wr>  -  (Nz/2)  (F|Vaa  +  2FaFbVab  +  F|Vbb) 
-  (Nz/2)  Fj(2VAB  -  VAA  -  VBB )s2 


Since  <Wr>  is  independent  of  temperature  in  this  approximation, 
by  equations  (10)  and  (12) 


<Wr>  -  E(s)  -  U(s) 


Thus ,  the  configurational  free  energy  becomes 
F(s)  -  F(0)  -  NKT[FA(FA+FBs)ln{FA(FA+FBs))  + 
fB(FB+FAs)ln{FB(FB+F*s)} 

+  2FAFB(l-s)ln{FAFB(l-s)}  -  2FAlnFA  -  2FBlnFB] 
-  NzVF| s2  . 


after  using  the  expression  for  g(s)  in  eq.  (5)  and  eq.  (15)  for 
the  energy  in  eq.  (11),  where  V  is  V'  as  defined  in  equation  1, 
in  the  nearest  neighbor  approximation.  For  the  alloy  of 
composition  FA  -  FB  -  1/2,  this  expression  reduces  to 

F(s)  -  F(l)  -  (NKT/2) ( (1+s) ln(l+s)  +  (l-s)ln(l-s)  -  21n2) 

+  l/4NzV ( 1 - s2 )  . (18) 


The  equilibrium  state  of  the  configuration  in  the  alloy  is  found 
by  minimizing  the  free  energy.  Thus,  the  condition 


gives 


s  -  tanh(zV/2KT)s 


This  equation  determines  the  equilibrium  value  of  s  as  a  function 
of  T.  The  above  equation  can  be  decomposed  into 


s  -  tanh  X  and  s  -  2KT/zNX  . 


When  these  two  equation  are  plotted  as  a  function  of  X  the 
intersection  of  the  curves  give  the  equilibrium  value  of  s  at  T 
(see  Figure  1).  There  are  two  solution  of  (21),  one  at  s-0  and 
one  at  s>0  for  all  T  less  than  some  T.  The  solution  s-0 
corresponds  to  the  maximum  free  energy  or  the  unstable  state  of 
order,  whereas  the  solution  of  s>0  corresponds  to  the  minium  free 


Determination  of  S 
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energy  of  equilibrium  value  of  order.  The  solution  of  s>0 
indicates  the  existence  of  the  ordered  state  at  temperature  below 
a  critical  value.  The  value  of  s  decreases  very  slowly  at  first, 
then  as  the  temperature  continue  to  increase  drop  off  abruptly  at 
Tc ,  and  vanishes  at  all  temperature  greater  than  Tc.  The 
critical  temperature  is  found  by  the  conditions 


—  s  —  0 


These  conditions  lead  to 


KTc  -  zV/2 


for  the  critical  temperature  of  an  alloy  with  composition  FA  -  F0 

-  1/2. 

A  higher  approximation,  developed  by  kirkwood  [29]  ,  is  one 
which  takes  into  account,  on  the  average,  the  eneregy  spread  of 
WT  for  a  given  s  around  its  avergy  value  <Wr>.  Although  it  is 
not  possible  to  know  the  actual  distribution  of  Wr  we  can  find 
the  moments  of  deviations  from  <Wr>.  We  write  equation  (10)  as 
U(s)/KT  -  ln[g(s)  J  -  ln[Z  e^r/KT]  . (24) 


and  expanding  the  exponentials , 


U(s)  <Wr> 

-  -  -lnl  -  -  +  ■ 

KT  KT 

g(s) 

where  <Wrn>  -t  (Wrn/g(s))  . 
r- 1 


1  <V2>  1  <dJ3> 

- - - +  ---(25) 

2!  (KT)2  3!  (KT)3 


Upon  further  expansion  of  the  logarithm  above,  U(s)  can  be 
written  as 


>  ^  .**  A  -  -  ^  •*  /•  / 
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the  specific  heats  do  arise.  The  critical  temperature  is  again 
found  by  the  condition  of  3F/3s  -  3zF/3s2  -  s  -  0,  and  is 


Kirkwood  solution  has  the  advantage  over  Bragg  and  Williams  in 
that  it  is  easily  adapted  to  different  types  of  lattices  and  to 
nonstoichiometry.  Also,  in  principle,  it  is  possible  to  obtain 
any  degree  of  accuracy;  however,  equation  (30)  tends  to  conerge 
very  slowly  at  temperatures  below  Tc  and  higher  order  terms  are 
necessary  for  reliable  results. 

Bethe's  solution  [28]  has  been  shown  to  be  Guggenheim's 
quasichemical  method  of  solution.  His  original  model,  as  will  be 
presented  here,  is  basically  a  construction  of  the  grand 
partition  function  for  a  small  group  of  lattice  sites  in  which 
the  variables  of  temperature,  volume  and  absolute  activities  play 
important  roles.  However,  the  original  model  was  not  constructed 
as  a  canonical  ensemble,  instead  Bethe  based  his  development  on 
nearest  neighbor  interactions  and  the  short  range  order 
parameter.  The  short  range  order  parameter,  as  defined  by  Bethe, 
is  not  based  on  a  and  y9  sites  in  the  lattice  but  on  the 
configuration  of  nearest  neighbors.  This  parameter  measures  the 
way  in  which,  on  the  average,  each  atom  is  surrounded  by  its 
neighbors,  or  the  extent  of  local  order.  The  short  range  order 
is  defined  as 

(q  -  qr) 

. . . (31) 

(qm  -  qr> 


.  4  *  %V 


\\" 

■  ’■  ‘-Ovn 


17 


where  q  -  Q^g/Q,  qm  is  the  maximum  value  of  q  and  qr  is  the  value 
of  q  in  the  disordered  state.  Here,  the  limits  of  zero  and  one 
are  obtained  for  complete  disorder  and  order  respectively,  as  for 
the  long  range  order. 

Now,  we  select  an  arbitary  group  of  sites  for  consideration 
on  the  basis  of  nearest  neighbor  interactions,  consider  for 
example  Figure  2.  The  group  consists  of  an  interior  and  a 
boundary,  all  remaining  sites  of  lattice  are  exterior.  The 
influence  of  the  exterior  atoms  are  taken  into  account  only  on 
the  average.  For  simplicity  the  method  is  developed  for  an  AB 
alloy. 

For  the  first  approximation  we  selected  an  a  site  as  the 
interior.  Thus  the  boundary  will  consist  of  z  nearest  neighbor  (3 
sites.  For  the  second  approximation  the  group  may  be  selected  as 
in  Figure  (2b) ;  the  larger  the  group  chosen,  the  greater  the 
accuracy.  However  the  mathematics  become  correspondingly  more 
complicated  [28].  Considering  the  first  approximations,  we 
assume  that  a  state  of  long  range  order  exists  in  the  exterior. 
This  ordering  will  affect  the  boundary  of  our  group  by 
encouraging  A  atoms  to  a  sites  and  B  atoms  to  d  sites.  This 
tendency  to  order  on  the  boundary,  due  to  the  exterior,  is 
denocci  by  an  ordering  energy  u.  The  energy  u  is  equal  to  the 
difference  in  energy  of  a  wrong  atom  on  an  a  or  f)  site  and  a 
right  atom  on  an  a  or  (3  site. 


Consider  now,  the  energy  of  the  alloy  which  is  determined  by 
the  number  of  pairs  QAA ,  QgB  and  QAB  for  a  definite  arrangement 
of  atoms  on  the  lattice  sites.  THe  energy  of  a  crystal  of  pure  A 
with  F*N  atoms  will  be  -VAA  (#  of  pairs)  or  -z/2  F*NVAA  -  - 
FaQVaa  and  similarly  for  a  crystal  of  B  atoms  only,  the  energy  is 
-FbVbbQ.  Thus  the  energy  of  the  alloy,  with  respect  to  the  pure 
crystals  of  its  components  is 

E  “  -VAA(feA  -  VabQab  -  VBBQ88  +  (VaaFaQ  +  VbbFbQ).  ----(32) 

This  reduces  to 

E  -  ( (1/2) (VAA  +  VBB)  -  7*810*8  -  *VQ*b  ----(33) 

since 

Qaa  “  EaQ  *  Qab/2 

Qbb  “  EBQ  *  Qab/2  •  (34) 

So  that  the  superlattice  will  be  stable  at  lower  temperature  V 
must  be  taken  less  than  zero.  The  Boltzmann  factor  is  then 

e-v/KT  „  e  .  (35) 

For  the  boundary,  the  relative  probablity  of  finding  an  A  atom 
instead  of  a  B  atom  due  to  the  exterior  is 

e-u/'T  -  x  .  (36) 

From  the  two  Boltzmann  factors  above  various  probabilities  in  the 
interior  and  in  the  boundary  may  now  be  found. 

If  the  central  atom  is  A,  a  right  atom,  the  relative 
probability  of  finding  n  wrong  A  atoms  in  the  boundary,  Pr(n) ,  is 

Pr(n)  -[^)«nxn  (37) 

If  the  center  atom  is  a  B  atom  the  probablity  of  finding  n  wrong 
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atoms  on  the  boundary  is 


Pw(n)  -  |nJxneZ"n  . (- 

since  each  of  the  z-n  right  atoms  in  the  boundary  has  an 
interaction  energy  V  with  the  wrong  interior  B  atom.  The 
binomial  coefficient  (2)  is  the  number  of  ways  of  arranging  n 
wrong  A  atoms  in  the  z  boundary  sites.  The  total  relative 
probability  for  the  interior  atoms  being  correct  is 


r,  “n5,Pr(n)  -  (1  +  ex)2 


and  for  being  wrong  is 


w,  -nZQPw(n)  -  ( €  +  x)2 
The  normalized  probabilities  are  then 


r~  -  - - -  ,  and  w_  -  -  .  . (41) 

a  r,+  w,  °  w^ 

It  is  also  necessary  to  caculate  the  relative  probability  of 

sites  in  the  boundary  being  wrongly  occupied.  This  is  the 

average  number  of  wrong  atoms  in  the  boundary  divided  by  z.  The 

relative  probability  of  finding  n  wrong  atoms  in  the  boundary  is 

Pr(n)  +  Pw(n) .  The  average  number  of  wrong  atoms  in  the  boundary 

is  calculated  as  follows: 


2„n(Pr(n)  +  Pw(n)) 

.  .  n—  0 

<n>  -  zwfl - i - 

Sn(Pr(n)  +  Pw(n) ) 

n 


---(42) 


Here,  w^  is  ‘.he  nomalized  probability  of  having  a  wrong  atom  in 
the  boundary  0  sites,  therefore 


x 


Now  in  order  to  determine  the  unkown  energy  u  in  x,  Be the 
introduced  a  symmetry  or  consistency  condition.  Since  there  are 
an  equivalent  number  of  a  and  £  sites  and  A  and  B  atoms  for  an  AB 
alloy  at  stoichiometry,  the  probability  must  be  symmetric,  i.e. 
the  same  probabilities  hold  whether  an  a  or  0  site  is  chosen  as 
the  origin  of  our  group .  Thus , 

wa  -  W£  . (44) 

This  is  called  the  consistency  equation.  Putting  equations  (41) 
and  (43)  into  the  above  equation  the  value  of  u  at  T  can  be  found 
as  a  function  of  V  and  T, 


x  +  e  1  1 

— : — -  -  x  .  . (45) 

1  +  ex 


The  short  range  order  a  and  the  configurational  energy  of 
the  alloy  can  now  be  found.  The  number  of  AB  pairs  is  z-n  for  a 
right  center  atom,  and  n  for  a  wrong  center  atom,  therefore  on 
the  average  the  fraction  of  AB  atoms  in  the  alloy  is 


Qa! 


q  - 


Z  Z 

(z-n)Pr (n)  +^SQnPw(n) 


-(46) 


Q  z(r,-  +  w,  ) 

For  the  AB  alloy  qffl  -  1  and  qr  -  0.5,  thus,  from  eq.  (31),  a  -  2q 
-  1,  or  from  equation  (46), 


,  4ex  1 

a  "  1  '  "  1+ex”  1  +  xZ/2-l 

The  configurational  energy  is  then 


(47) 


E  -  VQa8  -  -VQ,,  -  -  (VQ/2)  (a  +  1) 


(48) 


Measured  from  the  completely  ordered  alloy  this  energy  becomes 


E  -  E(0)  -  -  (VQ/2) (1  -  a). 


(49) 


The  long  range  order  can  be  found  by  use  of  equations  (2),  (39), 
(40) ,  and  (41)  and  is 


s  -  2ra  -  1 


(l+ex)z  -  (g+x)z 
(l+«x)z  +  («+x)z 


----(50) 


Recalling  that  at  complete  disorder 

ra  “  r/8  ”  wa  ”  w/9  “  1/2 

and  using  equation  (45) ,  the  critical  temperatures  for  an  AB 
alloy  is  given  by 


iSc  ln  A  |26'  '  . <5l> 

In  this  form,  extension  to  other  lattices  besides  simple  cubic  AB 
alloys  is  quite  laborious.  However,  this  is  equivalent  to  the 
quasichemical  method  of  Fowler  and  Guggenheim  [31];  and  Takagi's 
combinatorial  analysis  [32]  allows  a  simpler  extension  to 
different  compositions  and  structure. 

A  comparison  of  the  three  early  models  can  be  seen  in 
Figures  3,  4  and  5.  All  three  of  these  models  have  neglected 
important  physical  considerations,  as  do  the  more  recent  models, 
these  include:  (1)  interaction  of  atoms  which  are  not  nearest 
neighbors,  (2)  changes  in  cell  size,  lattice  symmetry  and  thermal 
expansion,  (3)  lattice  vibrations,  (4)  limit  of  the  validity  of 
the  ordering  energy  and  finally,  (5)  effect  of  size  differential 
in  the  ordering  phenomena.  A  comparison  of  theory  and 
experiment,  Figures  6  and  7  shows  poor  quantitative  agreement, 
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Figure  4.  Configurational  Energy  as  a  function  of  T.  [26] 


Figure  5.  Long-range  order  $  as  a  function  of  T.  [26] 


qualitatively  the  models  are  essentially  correct  [26],  The  more 
recent  models  are  somewhat  more  quantitatively  correct  and  will 
be  briefly  explained. 

Cluster-Variation,  Concentration  Waves  and 
Band  Theory 

The  cluster-variation  [39]  method  is  an  extension  of  the 
quasichemical  model.  This  model  uses  the  same  type  of  energy 
terms  as  the  earlier  theories.  However,  it  introduces  a  way  of 
arriving  at  a  better  approximation  to  the  entropy.  The  free 
energy  is  written  in  terms  of  the  different  configurations  a 
cluster  can  have  and  is  then  minimized  with  respect  to  the 
configurations.  For  example,  for  the  FCC  structure  a  tetragonal 
group  of  atoms  is  chosen  as  the  basic  cluster;  distribution 
variables  are  assigned  to  each  configuration  of  the  group  and  the 
free  energy  is  written  in  terms  of  these  variables.  Minimizing 
the  free  energy  with  respect  to  these  variables  produces  a  set  of 
simuitanious  algebraic  equations,  usually  of  higher  order.  The 
values  of  the  distribution  variables  which  satisfy  these 
equations  represent  the  equilibrium  configurations  of  the  system 
[40] .  The  early  formulation  of  the  method  did  not  allow  one  to 
choose  very  large  clusters  (which  increases  the  accuracy)  due  to 
the  large  number  and  nature  of  the  simuitanious  equations. 
Kikuchi,  however,  has  developed  an  iteration  scheme,  Natural 
Iteration  [39],  which  is  particularly  suited  to  the  method  and 
guarantees  convergence . 


A  number  of  authors,  [41,  42,  43],  have  used  this  method  and 
results  agree  fairly  will  with  experiment.  This  is  perhaps  one 
of  the  more  promising  methods  of  describing  cooperative  phenomena 

The  idea  of  concentration  waves  was  introduced  by  Landau 
[44]  and  Lifshitz  [45]  in  (1937,  1941)  in  the  development  of  the 
phenomenological  theory  of  phase  transitions  of  the  second  kind. 
This  method,  applied  to  the  order-disorder  transition  in  alloys 
[46]  has  proved  to  be  fruitful.  This  theory  enables  one  to 
understand  the  symmetry  aspects  of  the  order-disorder  transition, 
to  take  into  consideration  long  range  interactions  in  the  alloy 
and  to  predict  the  structure  of  the  ordered  phases,  which  the 
other  theories  do  not  allow,  based  on  concrete  interatomic 
potentials . 

All  of  the  preceeding  theories  reduce  to  the  comparison  of 
the  free  energy  of  the  ordered  and  disordered  phases  caculated 
from  nearest  or  at  best  next  nearest  interactions.  The  question 
of  whether  or  not  the  structure  of  the  ordered  phasee  is  stable 
[at  the  relevant  choice  of  interatomic  energies  i.e.  what  value 
of  i  is  chosen  in  equation  (1)]  and  what  the  structure  of  the 
ordered  phase  is  are  not  considered.  The  early  theories  as  well 
as  the  cluster- variation  inprovement  become  especially 
unsatisfactory  in  light  of  recent  developments  in  the 
pseudo -potential  theory  of  metal  alloys  and  covalent  compounds  (a 
large  part  of  the  C-Ti  bonds  are  covalent)  as  well  as  the  theory 
of  stress  induced  interactions  in  a  solid  solution.  Studies  show 


that  interatomic  interactions  which  are  responsible  for  ordering 
cannot  be  short  range.  Thus  for  substantial  improvement  of  the 
above  theories,  the  real  long  range  interactions  must  be  taken 
into  account. 

These  difficulties  are  overcome  by  using  the  method  of 
static  concentration  waves.  Atomic  distribution  in  a  binary 
alloy  can  be  described  by  means  of  the  function  n(t) .  This  one 
function  adequate  because  the  occupational  probabilities  nA (f ) 
and  ng(f)  for  A  and  B  atoms  respectively  are  not  indepedent. 

They  must  satisfy  the  identity 

nA(t)  +  ng(t)  -  1  ,  . (52) 

here  f  is  a  crystal  lattice  vector.  This  function  can  be  used 
for  interstitial  as  well  as  substitutional  alloys.  In  the 
disordered  state  the  probabilities  n(r)  are  the  same  for  all 
sites  which  can  be  occupied.  They  are  equal  to  the  atomic 
fraction  c  of  the  relevant  component  in  a  substitutional  alloy 
and  in  an  interstitial  solution  they  equal  the  fraction  of 
interstitial  sites  which  are  occupied.  In  the  ordered  phases 
n(f)  becomes  dependent  on  the  site  coordinate  T.  For  instance, 

if  n(t)  describes  t  values  n1  ,  n2  - nt  on  a  set  of  crystal 

sites  (f),  then  describes  the  sublattices  into  which  the 
disordered  lattice  has  been  subdivided  as  a  result  of  ordering. 
The  n(  are  occupational  probabilities  of  the  sites  of  the  1st, 

2nd  --  t-th  sublattice.  The  function  n(t)  can  be  expanded  in  the 
form  of  a  Fourier  series,  i.e.  it  can  be  represented  as  a 


superposition  of  static  concentration  waves: 

n(f )  -  c  +  1/2  2  [QCftjJexpCikj-r)  +  Q*(kj )exp( - ikj r) ] ,  (53) 

where  exp(ikr)  is  a  static  concentration  wave,  kj  is  a  non- zero 
wave  vector  defined  in  the  first  Brillouin  zone  of  the  disordered 
alloy,  t  is  a  site  vector  of  the  lattice  (f)  describing  a 
position  which  can  be  occupied  by  an  atom,  the  index  j  denotes 
the  wave  vector  in  the  Brillouin  zone  and  Q(kj )  is  a  static 
concentration  wave  amplitude.  Alternatively  this  can  be  written 
by  combining  those  terms  whose  wave  vector  kj's  enter  into  the 
stars  (the  star  is  a  set  of  wave  vector  kj  ,  which  may  be  obtained 
from  one  wave  vector  by  applying  to  it  all  operations  of  the 


symmetry  group  of  the  disordered  phase) : 

n(t)  -  c  +  2  t^EgU),  . (54) 

where  Es(r)  -  1/2  2  (r(js )exp( ifcj s r)  +  r*( js )exp( - ikj r)  (55) 

and 

Q(kj#)  -  n,rs(j9).  . (56) 


The  summation  is  carried  out  over  all  wave  vectors  in  the  star, 
n's  are  the  long  range  order  parameters  and  rs(js)  are 
coefficients  which  determine  the  symmetry  of  the  occupation 
probabilities  n(t)  with  respect  to  rotation  and  reflection 
symmetry  conditions.  The  requirements  that  in  the  completely 
ordered  state,  when  all  n(f)  are  either  zero  or  unity,  that  all 
n(t)  are  equal  to  unity  completely  defines  the  constants  r8(jg). 
It  turns  out  that  the  concentration  waves  are  related  to  the 
amplitude  of  coherent  scattering  in  X-Ray  diffraction  by  the 


3: 


equation 

Y(q)  -  1/2  (fB-fA)S  [Q(kj)S  exp(-i(q-kj )r  +  Q*(kj )Z  exp  - ( q+kj ) r ] 

1  . (57) 

where  fA  and  fB  are  the  form  factors  of  the  A  and  B  atoms 

respectively . 

Now,  the  configuration  part  of  the  Hamiltonian  of  the  alloy 


H  -  1/2  S  V(t,t')c(t)c(t'), 


where  the  summation  is  over  all  of  the  lattice  sites  and 
c(f)  -  1  if  there  is  a  solute  atom  at  r, 
c(r)  -  0  if  not. 

The  interchange  energy  V(r,t')  is  defined  by  the  equation 

V(r.r')  -vAA(f,t')  +vBB(t,t')  -  2Vab  (t , t ' ) .  (59) 

This  is  analogous  to  energy  defined  in  equation  (1),  however  it 
takes  into  account  all  interactions  in  the  lattice. 

The  problem  of  knowing  the  ordered  phase  structure  reduces 
to  determining  the  function  n(f).  Notice  that  n(f)  follows  a 
sort  of  Pauli  exclusion  principle:  each  lattice  site  is  either 
occupied  by  one  or  zero  atoms  of  some  definite  type.  Thus  we  can 
reason  that  it  follows  the  Fermi-Dirac  distribution: 


where  the  chemical  potential  n  is  determined  by  the  condition  of 
conservation  of  mass 


ri. 


where  N1  is  the  number  of  solute  atoms.  In  the  self-consistent 
field  approximation,  ^(f)  is 


*(t)  -  2  v<r.t')n(t). 


This  approximation  neglects  correlation  effects.  Substituting 
this  potential  into  (60) 


+  — S  V(r,f ' )n(f' )  +  1 


which  corresponds  to  the  free  energy 


fi  -  U  -  TS  -  n(t)  . (64) 

where  U  -  1/2  2  ,v(f ,r' )n(f )n(t* ) 

r  ,  r 

is  the  internal  energy  and 

S  -  -K2  [n(t)ln(n(£))  +  (l-n(t) )ln(l-n(£) ) ]  (65) 

is  the  entropy.  The  Helmholtz  free  energy  of  the  system  is 


F  -  U  -  TS. 


By  vising  the  expression  for  n(£)  in  equation  (54),  and  equation 
(63)  the  equilibrium  values  of  n,  at  T  can  be  found: 


c  +  2  n,E#(f)  -  exp 


-M  +  V(0)c  +  *2jv(fcs)n,E,(r) 


+  1.  (67 


Here  V(fc)  -  2  V(£)exp(ikr)  is  the  Fourier  transform  of  the 
interatomic  potential,  and  can  be  found  by  measuring  the  diffuse 
scattering  from  a  single  crystal  of  the  disordered  alloy,  [46] 


Idiff(k)  -  (fA-fB)2 


_ c(l-c) 

l  +  ---kt£'~ 


This  method  can  be  used  to  predict  the  critical  temperature, 
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ordered  structure  and  stability  apparently  with  much  greater  ease 
and  accuracy  than  the  former  theories  [46] .  The  critical 
temperature  calculated  for  the  CuAu3  alloy  by  this  method  is  only 
19  degrees  lower  than  the  actual  experimental  value.  This 
accuracy  is  exceptional  as  compared  to  the  previous  models .  It 
also  has  the  advantage  that  the  interaction  energies  are  in  no 
way  ambiguous  and  are  measurable  through  X-Ray  diffraction. 

One  of  the  band  theory  models  of  the  ordering  phenomena  [37] 
uses  a  tight  binding  Hartree  Hamiltonian  to  determine  the 
ordering  energy  incorporated  with  the  cluster -variation  method 
for  the  entropy  of  the  system.  Here  the  bulk  density  of  states 
is  derived  from  a  configurational  average  over  the  local 
densities  of  states  at  the  central  site  of  a  cluster.  The  free 
energy  of  the  alloy  is  determined  by  minimizing  a  model 
free-energy  function  over  the  space  of  possible  electron- ion 
configurations.  These  configurations  may  be  classified  by  the 
values  of  certain  order  parameters  [a  )  which  measure  multisite 
correlations,  e.g.,  long-range  order,  short  range  order, 
three-body  correlation,  etc.  Thus  at  a  temperature  T, 

F(T,c)  -  mina  F(T,c;a  ) 

-  E(T,c;a  )  -  TS(c,a  )  . (69) 

where  c  is  the  composition.  A  schematic  of  the  free  energy 
calculation  is  shown  in  Figure  8.  In  Figure  8,  Nina(w)  is  the 
local  density  of  states  and  <n,-n{7>  is  the  thermal  averaged 
occupation  of  sites. 
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Figure  8.  Block  Diagram  Showing  Basic  Steps  to  Self  Consistent 

Electronic  Theory  of  the  Alloys  Ordering  Energy  and  in 
Determing  Alloy  Thermodynamics. 


This  model  predicts  a  critical  temperature  of  the  CuAu 
system  which  is  in  error  by  19%  which  is  more  accurate  than  the 
early  set  of  models  but  can  not  compare  with  the  accuracy 
obtained  in  the  method  of  concentration  waves.  One  notable  point 
of  this  theory  however  is  that  it  actually  calculates  the  atomic 
interaction  erergies  as  a  function  of  short-range  order.  These 
energies  are  constant  only  near  a  -  0  and  can  vary  in  magnitude 
as  much  as  60%.  This  makes  the  early  Ising  models  and  the 
cluster-variation  model  even  more  unsatisfactory  since  they  are 
based  on  a  prior  knowledge  of  this  energy,  which  is  assumed 
constant  at  all  values  of  order. 

Since  Fosdick  [38]  first  applied  the  Monte  Carlo  method  to 
order-disorder  transformations  in  solids  a  number  of  attempts 
have  been  made  to  improve  and  extend  this  procedure  [47]. 
Unfortunately  the  stability  of  any  given  configuration  in  this 
model  is  again  based  on  the  nearest  neighbor  interaction 
energies.  Since  the  method  can  require  a  lot  of  computation  time 
and  exhibits  the  same  energetic  restrictions  as  early  Ising 
theories,  this  method  of  approaching  the  problem  of  ordering  does 
not  seem  to  be  worth  the  effort  in  this  research. 

Ordering  Energy 

It  seems  that  the  most  productive  model  available  to  predict 
the  ordering  phenomena  is  the  method  of  static  concentration 
waves.  This  is  the  theory  that  is  being  pursued  to  predict  the 


vv 


ordering  properties  of  the  C-Ti  system.  The  absence  of  ambiguity 
as  far  as  energy  terms  are  concerned  makes  this  theory  more 
desirable  than  the  Ising  model.  The  band  model  is  orders  of 
magnitude  more  involved,  as  well  as  less  accurate.  However,  it 
will  be  desirable  to  compare  predictions  made  by  the  Ising  models 
and  to  this  end  an  understanding  of  the  ordering  energy  in 
equation  (1)  is  necessary. 

Some  justification  of  choosing  nearest  neighbor  interaction 
energies  instead  of  long-range  effects  in  ordering  seems  to  be 
found  in  the  quantum  theory  which  holds  that  forces  between  atoms 
decrease  rapidly  as  the  distance  separating  them  increases.  As 
stated  earlier,  the  pair  interaction  energies  VAA ,  VBB  are  the 
interaction  energies  between  pairs  of  atoms  in  the  pure  elemental 
crystals  and  VAB  the  interaction  in  the  alloy.  However,  this 
definition  is  very  ambiguous.  What  do  we  measure  to  find  these 
energies,  or  how  do  we  calculate  them?  The  energy  V1  in  equation 
(1)  is  the  first  term  in  the  expansion  of  the  Fourier  Transform 
of  the  long-range  interaction  energy  used  in  the  static 
concentration  wave  model,  where  for  a  FCC  lattice,  [46] 

V(kg )  -  V(2*a*3)  -  -4w,  +  6w2  -  8w3  . (70) 


where  w  are  the  energies  corresponding  to  each  coordination  shell 
The  simple  nearest  neighbor  approximation  has  been  shown  to 
be  extremely  crude  [49,50],  for  instance,  it  cannot  explain  why 
most  order-disorder  transformations  are  first  order,  why  the 
maximum  value  of  Tc  often  occurs  away  from  stoichiometry  or  why 
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an  FCC  alloy  should  order  at  all.  Freidel  [50]  and  lthers  [51] 
have  been  able  to  show  that  because  conduction  electrons  in  a 
metal  are  limited  to  energies  below  the  Fermi  surface  a  complete 
screening  of  the  ionic  metal  cores  are  not  possible,  thus  long 
range  effects  are  present.  The  pseudo-potenial  theory  [52]  also 
points  to  many  body  interactions.  However,  neither  one  of  these 
methods  of  calculating  interaction  energies  can  be  used  for  the 
alloy  system  under  investigation  with  any  faith  due  to  the 
presence  of  covalent  bonding.  These  bonds  have  a  strong  angular 
nature  and  many  body  interactions  [23]. 

A  theory  recently  developed  to  find  the  ordering  energies  of 
transition  metal  alloys  [53]  is  available,  however  it  has  not 
been  used  extensively  and  requires  a  number  of  parameters  which 
might  be  as  difficult  as  ordering  energies  themselves  to  obtain. 

After  looking  in  several  references  [54,55,56]  we  found 
values  that  were  experimentally  verified  (see  table  1)  and  used 
these  in  calculating  the  critical  temperature  for  the  Cu-Au  alloy 
system.  The  critical  temperature  calculated  from  these  values 
using  the  Bragg-Williams  approximation  was  five  orders  of 
magnitude  higher  than  the  measured  value.  Clearly  the  heats  of 
formation  are  not  the  interaction  energies  which  we  are  looking 
for,  although  they  must  be  connected  in  some  way.  It  now  appears 
that  the  only  way  to  arrive  at  an  acceptable  number  for  these 
energies  is  through  X-Ray  diffraction. 
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Table  1 


Heats  of  Formation  of  Standard  States  from 
Monatomic  Gases 


Bond 

Energy 

Au-Au 

88.0  +  0.5  Kcal/mole 

Cu-Cu 

80.7  Kcal/mole 

Au-Cu 

55.4  t  2.2  Kcal/mole 

III.  EXPERIMENTAL  RESULTS  OF  PREVIOUS  INVESTIGATORS 


Phase  Diagram  and  Structure  of  C-Ti  Alloys 

Results  are  available  from  a  limited  number  of 
investigations  on  alloys  from  the  C-Ti  system.  The  phase  diagram 
for  the  C-Ti  system  is  shown  in  Figure  9.  Note  that  most  of  the 
lines  in  the  C-Ti  system  are  dotted.  This  indicates  the  lack  of 
adequate  data  to  firmly  establish  the  position  of  these  phase 
fields.  In  many  cases  these  dotted  phase  lines  are  based  upon 
calculation. 

Figures  10  through  16  [57]  show  the  microstructure  of  the 
C-Ti  alloys  with  different  carbon  contents  (0.4  w/o  to  17.5 
w/o) .  As  shown  in  Figure  13  to  16,  the  amount  of  7  in  the 
two-phase  region,  a  plus  7  (room  temperature)  ,  and  0  plus  7 
(temperature  over  920  degrees  C) ,  increase  with  increasing  carbon 
content  until  the  7  field  is  reached.  The  phase  within  the  7  in 
Figures  14  and  15  is  either  a  or  transformed  0  which  has 
precipitated  as  a  result  of  decreasing  solubility  in  TiC  with 
decreasing  temperature.  It  is  interesting  to  note  that  the 
surface  energy  of  the  7  in  contact  with  the  C-Ti  liquid  solution 
is  significantly  different  from  the  surface  energy  of  this  phase 
in  contact  with  the  0  phase  of  the  same  carbon  content.  In  the 
first  case  the  liquid  serves  as  the  matrix  with  7  as  the  included 
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C-Ti  Carbon  -Titanium 


Figure  9.  Phase  Diagram  for  the  C-Ti  Alloy  System. 


Figure io.  0.2  w/o  C,  6  hrs. 

at  920° C,  Water 
Quenched .  a , 
Transformed  0 
and  7 .  x200 


Figure  12.  11.2w/o  C,  1  min. 

at  1800° C,  Quenched. 
Liquid  Plus  7.  x200 


Figure  n.  0.2w/o  C,  as  Cast. 

7  Stringers  Plus 
Serrated  a.  x400 


Figure  13.  O.Uw/o  C,  36  hrs. 

at  800°C,  Water 
Quenched.  Equiaxed 
a  PLus  7-  x!50 


Figure  14.  4.0w/o  C,  48  hrs .  Figure  15. 
at  1200° C,  Water 
Quenched . 

Transformed  fi 
Plus  7.  Structure 
in  7  ts  8 


Figure  16.  17.5v/o  C.  120 
hrs.  at  900°C , 
Water  Quenched. 
Single  Phase  7. 


11 . 2w/o  C,  200 
hrs.  at  700° C, 
Water  Quenched 
7  Plus  a. 
Structuer  in  7 
is  a.  x500 


phase;  and  in  the  second  case  the  7  tends  to  serve  as  the  matrix 


for  the  as  is  indicated  by  the  precipitates  in  the  7  phase 

and  by  the  sharp  angles  the  7  constituents  have  developed  at 
grain  boundaries  of  f3 ,  Figure  14.  Figure  15  gives  the  appearance 
of  the  liquid  serving  as  matrix  around  the  phase  for  a  different 
composition  alloy.  The  7  phase  has  a  face-centered  cubic  X-ray 
diffraction  pattern  and  the  compound  was  a  NaCl  type. 

Structure  and  Properties  of  Low  Carbon  Weight 
Percentage  C-Ti  alloys  (C<0.5w/o) 

Previous  investigations  in  the  C-Ti  system  have  concentrated 
on  determining  and  understanding  the  properties  of  low  weight 
percentage  C-Ti  alloys  (C<0.5w/o)  [58-59],  single  crystal  TiC 
[60-63]  and  ploycrystalline  TiC  [64-68].  The  properties  of  the 
three  binary  C-Ti  alloys,  given  in  Table  2  [59],  illustrate  the 
solid-solution  strengthening  effect  of  carbon.  Figure  17  shows 
the  effect  of  heat-treating  temperature  on  the  properties  of  a 
Ti-0.16w/o  C  alloys.  Microstructures  for  corresponding 
conditions  are  given  in  Figure  18 .  When  this  alloy  is  annealed 
and  quenched  from  the  /9-carbide  field,  there  is  a  marked  decrease 
in  strength  and  toughness  caused  by  the  formation  of  TiC,  and 
consequent  lowering  of  the  carbon  in  solid  solution.  Data  for 
the  Ti-0.27w/o  C,  the  strength,  particularly  the  yield  strength, 
is  increased  by  dissolving  the  TiC  particles  through  the  use  of 
high  temperature  a  anneal.  The  high  tensile  ductilities  are 
found  when  the  Ti-0.27w/o  C  alloy  is  quenched  from  the  a-/3 


Properties 


traiufonned  to  a  on  qucnchln*. 
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field.  At  a  level  of  0.47w/o  C  the  strength  properties  are 
greatly  increased  when  the  carbon  is  put  in  solid  solution 
through  the  use  of  increasingly  higher  annealing  temperatures, 
and  they  are  decreased  sharply  when  the  alloy  is  quenched  from 
the  ^-carbide  field.  Figure  21  shows  the  properties  of  a 
Ti-0.47w/o  C  alloy.  Microstructures  of  this  alloy  are  given  in 
Figure  22. 

A  further  illustration  of  how  heat  treatment  can  affect  the 
properties  of  Ti-C  alloys  is  given  in  Figure  23,  where  properties 
are  plotted  vs  carbon  content  for  three  different  heat 
treatments.  Tensile  ductilities  and  toughness  decrease  with 
increasing  carbon  content,  whether  in  solution  or  not.  Carbon  in 
solid  solution  appears  to  decrease  tensile  ductility  and  increase 
toughness  slightly. 


Mechanical  Properties  of  Titanium  Carbide 
Polycrystalline  TiC  have  been  produced  by  hot  press 
sintering  of  near  stoichiometric  TiC  powders.  The  results  show 
that  TiC  is  brittle  at  room  temperature  and  becomes  ductile  above 
800  degrees  C.  Katz  and  Lipsitt  [64]  have  shown  these 
stoichiometric  compositions  of  TiC  to  be  highly  sensitive  to 
grain  boundary  impurities.  Above  1000  degrees  C  a  grain  boundary 
segregant  containing  Si,  Fe,  Cu  becomes  near  molten  and  results 
in  greatly  reduced  strength.  A  large  amount  of  porosity  has  been 
identifies  in  scanning  electron  micrographs  of  the  fractured 
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Figure  21.  Properties  of  a  Ti-0.47w/o  C 

Alloy,[59] 
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Figure  22.  Microstructures  of  a  Ti-0.47v/o  C  Alloy. 
x250 [59] 


specimens  and  is  attributed  to  the  heat  treatment  and  fabrication 
process.  This  porosity  is  expected  to  adversely  affect  the 
material  ductility  and  strength. 

The  optical  micrographs  of  polycrystalline  substoichiometric 
TiC  with  different  compositions  are  shown  in  Figure  24  [67], 

Grain  sizes  ranged  from  14  to  22  jum  and  increased  monotonically 
with  decreasing  stoichiometry.  Pores  were  evenly  distributed, 
and  also  grew  with  decreasing  stoichiometry.  Hardness  was  shown 
to  drop  steadily  as  expected  with  decreasing  C/Ti  atom  ratios 
(Figure  25) . 

Compression  Test  Results 

Typical  nominal  stress -nominal  strain  curves  (o-O  [68]  for 
1700  degrees  C  annealed  polycrystalline  TiC,  having  an  average 
grain  size  of  8.2  pm,  deformed  in  compression  at  900  degrees  C  to 
1700  degrees  C  are  shown  in  Figure  26.  Yield-point  behavior  is 
observed  for  specimen  deformed  at  1000  degrees  C.  The  yield 
stress  in  polycrystalline  TiC  is  strongly  dependent  on  the 
temperature  of  deformation,  as  shown  in  Figure  27.  Similar 
observations  were  reported  for  both  single-crystal  [60]  and 
polycrystalline  TiC  [64,67],  The  dependence  of  the  yield  stress 
on  grain  size,  as  shown  in  Figure  28,  follows  a  Hall-Petch  type 
[68]  relation 

<7y  -  o0  +  kd-’/2  . (71) 

where  ay  is  the  yield  stress,  ct0  a  friction  stress,  K  a  constant 
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Figure  27.  Yield  Stress  as  a  Function  of 
Temperature  in  Polycrystalline 
TiC  Deformed  in  Compression.!  68] 


Figure  28.  Yield  Stress  vs  (Grain  Size)*1^ 
for  Polycrystalline  TiC 
Deformed  in  Compression  at 
Temperature  Indicated.[68] 
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associated  with  the  propagation  of  deformation  across  the  grain 
boundaries,  and  d  the  grain  diameter. 

Figure  29a  [66]  is  an  SEM  micrograph  of  the  fracture  surface 
of  a  specimen  tested  at  room  temperature.  The  fracture  is 
predominantly  transgranular  cleavage.  Figure  29b  is  the 
corresponding  fractograph  for  a  specimen  tested  at  1100  degrees 
C.  Here  the  fracture  is  much  nore  intergranular,  and  the  extent 
of  intergranular  fracture  continues  to  increase  as  the  test 
temperature  is  raised.  Similar  observation  was  reported  for 
polycrystalline  TiC  [68]. 


Deformation  Mechanisms 

In  Tic,  three  thermally  activated  dislocation  mechanisms 
[68]  have  been  proposed  to  operate  during  deformation  at 
different  test  temperatures:  (1)  Peierls  stress  [69]  (<1200°C); 
(2)  diffusion  of  carbon  [64]  (1200°C  to  1400°C);  and  (3) 
diffusion  of  titanium  [65]  (>1500°C).  A  phenomenological 
relation  commonly  used  to  establish  a  law  governing  the 
deformation  behavior  is  the  dependence  of  the  yield  stress  (cry) 
on  temperature,  which  can  expressed  as 

<7  -  <70 exp (B/KT)  . (72; 


where  is  a  constant,  K  the  Boltzman  constant,  B  an  activation 
energy,  and  T  the  absolute  temperature.  Figure  30  is  a  plot  of 
logoy  as  a  function  of  1/T  for  polycrystalline  TiC  having 
different  grain  sizes  tested  in  compression.  A  change  in  slope 


Figure  29.  SEM  Fractograph  of  TiC  Specimens  Hot  Pressed  from 
Unsieved  Starck  Powder.  Bending  Tests  Performed  at 
(a)  Room  Temperature,  (b)  1100°C.[68] 
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Figure30.  Log  of  Yield  Stress  vs  Reciprocal 
of  Absolute  Temperature  for 
Polycrystalline  TiC  Deformed 
in  Compression. [ 68 ] 


was  observed  at  1500  degrees  C.  Below  1500  degrees  C  the 
mangnituta  of  the  slope  is  nearly  the  same  for  specimens  having 
four  grain  sizes;  above  1500  degrees  C  the  slope  is  larger  by 
more  than  a  factor  of  3.  Other  investigations  [60,70]  in  both 
single-crystal  and  polycrystalline  TiC  have  indicated  the 
existence  of  one  such  change  in  slope.  However,  the  temperature 
at  which  this  change  occurred  was  determined  to  be  different  by 
various  investigators.  It  has  also  been  found  to  depend  on  the 
degree  of  substoichiometry  of  TiC  [60], 


IV.  EXPERIMENTAL  PROCEDURES 


that  contains  a  tantalum  sheet  encircling  the  inside  diameter. 
The  crucible  was  supported  by  a  ceramic  rod.  The  induction  coil 
encircles  the  outside  of  the  quartz  tube.  The  tantalum  sheet 
serves  as  a  radiation  reflector  but  additionally  serves  as  a 
gathering  material  to  protect  the  melt  from  oxidation.  The 
graphite  crucible  is  the  susceptor  and  quartz  tube  serves  as  the 
containment  of  the  inert  gas  atmosphere.  Both  the  graphite 
crucible  and  graphite  powder  were  spectrographic  grade  material. 
The  bottom  ceramic  rod  was  used  to  limit  heat  transfer  by 
conduction  out  of  crucible.  An  infrared  sensor  was  used  to 


monitor  the  temperature  of  the  melt.  This  sensor  is  suspended 
from  the  top  support  of  the  heating  system.  Figures  31  and  32 
show  the  system  in  operation.  An  argon  atomosphere  was 
maintained  in  the  quartz  tube  during  heating  and  cooling  of  the 
melts.  A  graphite  bar  was  used  to  stir  the  fluid  alloy  mixture 
to  increase  the  diffusion  of  carbon  into  titanium,  enhancing  the 
formation  of  a  homogeneous  material  composition  and  reducing 
porosity  of  the  melts .  The  alloys  in  the  0  to  7  w/o  carbon 
region  were  produced  using  an  induction  generator  set  to  provide 
a  direct  current  of  40  amps  at  230  volts.  The  melts  were  limited 
in  size  to  0.35  inches  in  diameter  and  0.5  inches  high. 

Alloys  with  compositions  greater  than  7  w/o  carbon  were 
produced  by  arc  melting  under  an  inert  gas  atomosphere.  A  gas 
cooled  electrode  which  is  composed  of  Tungsten- 2%  Thoria  was  used 
to  produce  the  melts.  Initially,  problems  were  incountered  when 
the  cooling  gas  blew  the  graphite  powder  out  of  the  crucible.  A 
method  was  developed  in  which  the  graphite  powder  was  sealed  into 
holes  that  were  drilled  in  the  pure  Ti  rod.  This  enabled  high 
carbon  weight  percent  alloys  to  be  prepared.  Quality  of  alloys 
produced  by  the  arc  melting  process  was  less  desirable  than 
alloys  produced  by  induction  melting.  Arc  melted  alloys  tended 
to  be  less  homogenous  and  contained  amount  of  porosity. 


Fabrication 


Techniques  have  been  developed  for  the  fabrication  of 
subsize  compression  specimens  [71]  from  C-Ti  alloys.  Because  of 
the  small  size  and  properties  of  the  available  material,  a 
diamond  impregnated  soft  core  tube  (Figure  33)  which  has  an 
outside  diameter  of  6.4  mm  and  inside  diameter  of  4.4  mm  was  used 
to  prepare  cylindrical  compression  specimens.  The  drill  utilizes 
water  as  a  lubricant  and  coolant  to  protect  the  samples  from 
overheating.  After  a  cylinder  is  cut  from  the  melt,  a  diamond 
slitting  saw  was  used  to  cut  both  ends  of  the  cylinder  to  length 
and  obtain  parallel  ends.  Low  carbon  content  alloys  were  gummy 
in  nature,  thus  requiring  sharpening  of  the  diamond  blade  with  a 
carbon  block  to  rmove  cut  particles  stuck  in  the  blade.  Figure 
34  shows  the  diamond  slitting  saw.  Although  the  alloys  were  not 
excessively  hard,  large  amounts  of  time  were  required  to  prepare 
the  cylindrical  compression  sample,  which  measured  4.40  mm  in 
diameter  by  6.86  mm  high.  A  diameter  to  length  ratio  of  1.50  to 
1.75  was  used  for  all  tested  specimens.  These  specifications 
meet  ASM  Metal  Handbook  requirements  for  subsize  compression 
samples  [71],  Prior  to  compression  testing,  specimens  were 
ultrasonically  cleaned  in  acetone.  Figure  35  shows  a  typical 
compression  specimen. 
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Fabrication 


Techniques  have  been  developed  for  the  fabrication  of 
subsize  compression  specimens  [71]  from  C-Ti  alloys.  Because  of 
the  small  size  and  properties  of  the  available  material,  a 
diamond  impregnated  soft  core  tube  (Figure  33)  which  has  an 
outside  diameter  of  6.4  mm  and  inside  diameter  of  4.4  mm  was  used 
to  prepare  cylindrical  compression  specimens.  The  drill  utilizes 
water  as  a  lubricant  and  coolant  to  protect  the  samples  from 
overheating.  After  a  cylinder  is  cut  from  the  melt,  a  diamond 
slitting  saw  was  used  to  cut  both  ends  of  the  cylinder  to  length 
and  obtain  parallel  ends.  Low  carbon  content  alloys  were  gummy 
in  nature,  thus  requiring  sharpening  of  the  diamond  blade  with  a 
carbon  block  to  rmove  cut  particles  stuck  in  the  blade.  Figure 
34  shows  the  diamond  slitting  saw.  Although  the  alloys  were  not 
excessively  hard,  large  amounts  of  time  were  required  to  prepare 
the  cylindrical  compression  sample,  which  measured  4.40  mm  in 
diameter  by  6.86  mm  high.  A  diameter  to  length  ratio  of  1.50  to 
1.75  was  used  for  all  tested  specimens.  These  specifications 
meet  ASM  Metal  Handbook  requirements  for  subsize  compression 
samples  [71] .  Prior  to  compression  testing,  specimens  were 
ultrasonically  cleaned  in  acetone.  Figure  35  shows  a  typical 
compression  specimen. 


and  the  maximum  output  voltage  of  both  load  and  stroke  cells  is 
10  volts,  Calibrated  subscales  of  10%,  20%,  and  50%  are 
available  for  both  stroke  and  load. 

Data  Acquisition  Unit 

The  data  acquisition  system  (Figure  36)  is  composed  of  three 
parts  :  an  HP-7046B  X-Y  recorder,  a  Nicolet-3091  digital 
oscilloscope  with  bubble  memory,  and  an  HP-3497A  data  acquisition 
unit  combined  with  an  HP- 9836  minicomputer.  These  data 
acquisition  systems  monitor  the  compressive  data  coming  from  the 
load  and  stroke  sensors,  and  transmit  the  data  to  the  computer 
which  analyzes,  plots,  and  prints  the  incoming  data. 

Heating  Elements 

Compression  tests  were  conducted  at  elevated  temperatures 
(159  degrees  C  -  850  degrees  C)  using  a  specially  fabricated 
fumance .  The  furnace  consisted  of  a  quartz  tube  for  electrical 
insulation  surrounded  by  30  turns  of  nichrome  wire  1.0  mm  in 
diameter.  The  wire  was,  in  turn,  covered  by  a  layer  of 
fiberglass  which  served  as  a  thermal  insulator  (Figure  37) .  The 
dimensions  of  the  quartz  tube  were  30*32*160  mm.  Power  to  the 
heating  element  was  supplied  by  a  variac  transformer,  with  the 
temperature  controlled  by  adjusting  the  output  voltage  of  the 


transformer . 


Experimental  Procedure 

Figure  38  shows  the  flow  chart  of  the  MTS  data  acquisition 
system.  Chosen  strain  rate  was  obtained  by  setting  the  maximum 
stroke  length  and  adjusting  the  time  to  reach  the  total  stroke 
length.  The  relationship  is  shown  below  : 

Maximum  Ls 

strain  rate  -  - 

(Lg)  *  (Time  To  Reach  Maximum  Ls) 


Room  Temperature  Compression  Tests 
The  strain  rate  was  fixed  at  10*4  sec*1  for  room  temperature 
compression  tests.  The  X-Y  recorder  and  HP-3497A  data  acqusition 
unit  were  used  seperately  to  continuously  digitize  and  store  the 
data  from  the  load  and  stroke  sensors  of  the  MRS  machine.  For 
these  tests,  the  range  of  load,  coarse  stroke,  and  fine  stroke 
were  set  at  100  KN,  10  mm,  and  100%  respectively.  Material  type, 
specimen  number,  gage  length,  cross-sectional  area  and  specifics 
of  each  test  specimen  were  entered  into  the  computer  before 
performing  the  test.  The  specimen  was  placed  on  the  compression 
test  platens ,  which  were  attached  to  the  ram  of  the  MTS  machine . 
Upon  finishing  the  tests,  the  computer  plotted,  printed  and 
stored  the  data  automatically.  The  X-Y  recorder  was  used  to  plot 
the  stress-strain  diagram  while  the  test  was  in  progress. 


MTS  Test  Machine 


Flow  Chart  of  MTS  Data  Acquisition 


Elevaced  Temperature  Compression  Tests 


The  specimen  was  placed  between  cylindrical  linkages  made  of 
Inconel  706,  and  the  quartz  tube  wound  heating  element  placed 
over  the  linkages.  This  assembly  was  then  attached  to  the  load 
frame.  The  furnace,  supported  by  a  three-prong  clamp  (Figure 
39),  was  flushed  continously  with  argon  during  testing  to  prevent 
oxidation  of  the  specimen.  Before  testing,  the  specimens  were 
heated  for  20  minutes  to  allow  thermal  equilibrium  to  be 
reached.  The  chosen  temperatures  (150°C,  350° C,  450° C,  550° C, 
and  850° C)  were  controlled  and  maintained  by  varing  the  output 
voltage  of  a  variac  transformer.  Temperature  was  measured  using 
K-type  thermocouples  mechanically  attached  to  the  central  gage 
section  of  the  specimen. 

Failure  Analysis 

To  investigate  the  deformation  and  fracture  mechanisms, 
selected  specimens  were  examined  in  detail  with  an  JSM-840 
scanning  electron  microscope  operated  at  20KV.  The  scope  is 
equipped  with  Energy  Dispersive  Spectrometers  (EDS)  capable  of 
compositional  analysis.  The  SEM  examination  involved  the  study 
of  both  macroscopic  and  microscopic  surface  features. 

Macroscopic  features  were  studied  from  20X  to  50X.  The  center 
and  edge  of  the  fracture  surface  were  examined  for  features.  All 
microscopic  features  were  studied  from  200X  to  600X. 

Additionally,  the  microstructure  of  the  C-Ti  alloys  was  examined 
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using  a  Zeiss  photomicroscope.  Specimens  for  optical  microscopy 
were  polished  and  etched  with  a  mixture  of  water,  Nitric  Acid, 
and  Hydrofluoric  Acid  in  proportions  of  5:5:2  for  1  second. 
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v.  EXPERIMENTAL  RESULTS 


Optical  Metallography 

Metallographic  samples  were  prepared  from  the  bottom  surface 
of  the  compression  sample.  Figures  40  though  45  display 
representative  micro structures  of  the  C-Ti  alloys.  Samples  were 
produced  by  the  previously  mentioned  induction  melting  procedure 
or  arc  melting  procedure.  Figure  40  is  the  photomicrograph  of  a 
2  weight  percent  (w/o)  C-Ti  alloy.  Figures  41  through  45  are 
microstructures  of  alloys  with  increasing  carbon  contents  (4w/o 
to  lOw/o) .  The  microstructures  consist  principally  of  a  mixture 
of  Alpha-Ti  and  Gamma  phases  at  room  temperature.  A  third  carbon 
rich  phase  was  found  in  some  of  the  higher  carbon  percentage 
alloys  and  is  thought  to  be  the  result  of  localized  incomplete 
diffusion  of  carbon  into  the  melt.  A  comparison  of  the  results 
with  the  published  phase  diagram  for  the  C-Ti  system  (Figure  9) 
show  agreement  with  the  anticipated  phases  and  the  microstructure 
for  the  alloys  of  2  to  10  w/o  carbon  with  the  above  noted 
exception  of  the  high  carbon  content  phase . 

Figures  46  though  51  show  the  Gamma-phase  particle  size 
distribution  and  the  volume  fraction  of  Gamma  phase  (2-10  w/o 
C) .  The  Gamma-phase  particle  size  distribution  was  calculated  by 
using  the  intercept  method.  Five  microphotographs  of  the  same 


Figure  46.  Gamma  Phase  Size  Distribution  :  Ti-2w/o  C.  Microstructure  of 
Alloy  shown  In  Figure  40. 
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Figure  50.  Gamma  Phase  Size  Distribution  :  Tl-7w/o  C.  Microstructure  of 
Alloy  shown  in  Figure  44. 
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specimen  were  taken  prior  to  the  compression  test  to  caculate  the 
volume  fraction  of  Gamma -phase  of  the  alloys.  The  phase 
distribution  was  measured  using  a  HP-9111A  graphics  tablet  in 
conjunction  with  an  HP-9836  microcomputer.  The  mean  value  of 
five  measurements  from  each  specimen  was  included  in  this  study. 
The  results  show  that  the  volume  fraction  of  Gamma  phase 
increases  as  the  carbon  content  is  increased.  Because  of  this, 
the  strength  of  C-Ti  alloys  increases. 


Hardness  Measurements 

Hardness  measurements  were  made  on  the  alloys  to  aid 
characterization  of  their  mechanical  properties.  The  macroscopic 
hardness  measurements  were  found  to  vary  with  carbon  content. 

The  rockwell  C  hardness  varied  from  a  value  of  22  for  a  2  w/o 
carbon  alloy  to  40  for  the  10  w/o  carbon  alloy.  In  addition  to 
macrohardness  measurement,  microhardness  measurements  (using  a 
Vicker  hardness  unit)  were  made  on  the  Alpha-Ti  and  Gamma  phase 
structures  seen  in  the  photomicrographs  of  figures  40  though  45. 
The  diamond  indentations  on  the  Alpha-Ti  phase  were  clear  and 
large,  but  due  to  the  increased  hardness  of  Gamma  phase 


j* 


Si 


indentations  appear  much  smaller.  Figure  52  and  figure  53 
compare  the  microhardness  of  Alpha-Ti  and  Gamma  phase  for  the  2 
w/o  C  and  10  w/o  C  alloys.  From  this  data  it  is  concluded  that 
the  microhardness  of  the  Gamma  phase  increases  as  the  carbon 
content  is  increased.  The  Vicker  hardness  value  of  Alpha-Ti  is 


85 


insensitive  to  carbon  content.  The  hardness  of  C-Ti  is  therefore 
dictated  by  the  volume  fraction  and  distribution  of  Gamma  phase. 

Compression  Test  Results 
Room  Temperature  Compression  Tests 
Compression  tests  at  a  fixed  strain  rate  (10'4  sec'1)  were 
performed  on  the  C-Ti  alloys  at  differing  temperatures,  with  a 
total  of  34  specimens  tested.  Tabulated  results  of  all  room 
temperature  compression  tests  are  given  in  APPENDIX  C,  Figure  54 
shows  a  typical  compression  specimen  after  testing,  Macroscopic 
failure  occurred  along  a  45  degree  plane.  The  compression  test 
results  will  be  discussed  in  the  following  section,  Yield  Stress 
and  True  Failure  Strain. 

Yield  Stress  (Y.S.) 

A  summary  of  the  mechanical  properties  of  some  of  the 
tested  C-Ti  alloys  at  room  temperature  is  given  in  Table  3 . 

Figure  55  shows  the  characteristic  stress  vs  strain  curve  for 
C-Ti  alloys  with  different  volume  fraction  of  Gamma  phase.  This 
curve  indicates  that  the  elastic  modulus  of  C-Ti  alloys 
increases  as  the  volume  fraction  of  Gamma  phase  increases  Th< 
plot  of  yield  stress  as  a  function  of  Gamma  phase  fraction 
shown  in  Figure  56.  The  data  points  for  100%  Gamma  ;  •  • 

are  taken  from  references  [67]  and  [  T  '  Yield  $•:*■■■•  .. 

obtained  by  using  the  0  2%  offset  v :  e  1  d  .  • 
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Table  3 

Summary  of  Average  Mechanical  Properties 
for  C-Ti  Alloys  at  Room  Temperatuer 


Specimen 

No. 

Y.S.  True  M.F.S. 

(Ksi)  (Ksi) 

True  Failure 
Strain  (%) 

%Gamma 

Phase 

W1 

166 

189 

1.8 

80% 

118 

103 

166 

7.8 

60% 

17 

90 

176 

14.7 

40% 

128 

60 

154 

22.3 

20% 

127 

50 

143 

26.0 

10% 

Ti* 

20 

34 

54.0 

0% 

* - 

Data  is  taken 

from  Ref.  [72] . 

True  Strain  (in/in) 

Characteristic  True  Stress  vs  True  Strain  Curve  for  C-Ti 
Alloys  with  Different  Volume  Fraction  of  Gamma  Phase. 
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Figure  56.  Yield  Stress  as  a  Function  of  Gamma  Phase, 


points  indicate  a  dramatic  rise  exponentially  in  yield  stress  as 
the  volume  fraction  of  Gamma  phase  increases.  The  values  ranged 
from  168  Ksi  for  82%  Gamma  phase  to  58  Ksi  for  10%  Gamma  phase. 


True  Failure  Strain  (£f) 

True  failure  strain  was  measured  by  the  following  method. 
The  true  failure  plastic  strain  on  the  true  stress-true  strain 
curve  was  caculated  by  drawing  a  line,  with  slope  equal  to  the 
elastic  modulus,  from  the  failure  point  on  the  curve  to  the 
intersection  with  the  strain  axis.  This  intersection  denotes 
true  failure  plastic  strain  of  the  specimen.  A  plot  of  true 
failure  plastic  strain  as  a  function  of  the  volume  fraction  of 
Gamma  phase  is  shown  in  Figure  57.  There  is  a  dramatic  decrease 
in  the  ductility  as  the  volulme  fraction  of  Gamma  phase  rises. 
The  results  are  very  promising  in  that  alloys  containing  82%  and 
95%  Gamma  phase  still  show  2-5%  ductility.  This  indicates  that 
there  must  be  a  significant  fraction  of  the  Gamma  phase  which  is 
deforming. 


Elevated  Temperature  Compression  Tests 
Appendix  D  lists  the  mechanical  properties  of  the  C-Ti 
alloys  which  were  tested  in  compression  at  temperatures  from  150 
to  850  degrees  C.  A  total  of  thirty-four  specimens  were  tested 
in  the  manner  discribed  in  Section  IV.  The  results  are  presented 
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in  the  following  sections,  True  Maximum  Flow  Stress,  Yield 
Stress,  and  True  Failure  Strain. 


True  Maximum  Flow  Stress  (M.F.S.) 

The  true  maximum  flow  stress  is  directly  obtained  from  a 
true  stress -true  strain  curve.  Figures  58  and  59  show  the 
characteristic  true  stress  vs  true  strain  curve  for  C-Ti  alloys 
in  different  temperatures  ranging  from  150  to  850  degrees  C. 
These  curves  indicate  that  the  maximum  flow  stress  decreases  as 
the  test  temperature  increases  for  a  constant  volume  fraction  of 
Gamma  phase .  A  summary  of  the  mechanical  properties  of  chosen 
tested  alloys  for  57%  &  10%  Gamma  phase  at  elevated  temperatures 
is  given  in  Table  4.  A  plot  of  true  maximum  flow  stress  as  a 
function  of  Gamma  phase  is  shown  in  Figure  60.  The  true  maximum 
flow  stress  increases  exponentially  as  the  volume  fraction  of 
Gamma  phase  increases  for  a  constant  test  temperature. 


True  Failure  Strain  and 
Yield  Stress 

The  method  of  measuring  yield  stress  and  true  failure  strain 
has  been  described  previously.  Figure  61  shows  the  plot  of  yield 
stress  as  a  function  of  Gamma  phase.  The  yield  stress  increases 
as  the  test  temperature  decreases  and  drops  exponentially  with 
decreasing  volume  fraction  of  Gamma  phase  for  a  constant 
temperature.  The  dependence  of  yield  stress  as  a  function  of 


True  Strain  (In/In) 

e58.  Characteristic  True  Stress  vs  True  Strain  Curve  at  Elevated 
Temperatures  for  57%  Gamma  Phase. 


Table  4 
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Yield  Stress  as  a  Function  of  Volume  Fraction  of  Gamma  Phase 
at  Elevated  Temperature. 


fraction  of  Gamma  phase  nearly  are  identical  at  different  test 
temperatures . 

The  plot  of  true  failure  strain  as  a  function  of  fraction  of 
Gamma  phase  is  shown  in  Figure  62.  The  true  failure  plastic 
strain  decreases  exponentially  with  increasing  volume  fraction  of 
Gamma  phase  for  constant  temperature.  The  data  points  indicate 
that  the  slope  becomes  larger  as  the  test  temperature  falls  below 
350  degrees  C. 

Related  The  Results  At  Room  Temperature 
and  Elevated  Temperatures 

Figures  63  through  65  display  the  true  maximum  flow  stress, 
yield  stress,  and  true  failure  strain  as  a  function  of  reciprocal 
test  temperature  with  different  volume  fraction  of  Gamma  phase. 
The  data  points  for  the  100%  Gamma  phase  and  Ti  are  taken  from 
references  [67],  [68],  [72],  and  [73].  A  change  in  the 
controlling  deformation  mechanism  is  seen  to  occur  as  indicated 
by  the  change  in  slope  of  the  plot.  The  changes  in  slope  occur 
at  different  transition  temperatures  exist  for  C-Ti  alloys 
containing  differrent  volume  fraction  of  Gamma  phase.  This 
indicates  that  the  changes  in  controlling  mechanism  is  a  function 
of  both  temperature  and  composition.  The  range  of  transition 
temperature  is  from  470  degrees  C  for  10%  Gamma  phase  to  500 
degrees  C  for  60%  Gamma  phase.  Figure  66  plots  the  transition 
temperature  as  a  function  of  volume  fraction  of  Gamma  phase .  The 
transition  temperature  increases  linearly  with  increasing  volume 


Maximum  Flow  Stress  as  a  Function  of  Temperature. 


Figure  65.  True  Failure  Strain  as  a  Function  of  Temperature. 
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Figure  §6  Transition  Temperature  as  a  Function  of  Volume  Fraction  of  Gamma  Phase 


transition  temperature  increases  linearly  with  increasing  volume 
fraction  of  Gamma  phase. 


Scanning  Electron  Microscopy 
To  investigate  the  macroscopic  and  microscopic  fracture 
surface,  two  sets  of  specimens  tested  at  room  temperature  and 
elevated  temperatures  were  selected  for  SEM  study.  Each  set 
contained  four  specimens  having  the  same  volume  fraction  of  Gamma 
phase,  57%  and  10%.  All  specimens  were  cleaned  in  an  Acetone 
solution  to  remove  the  oxide  layer  on  the  fracture  surface.  The 
examination  procedure  is  described  in  detail  in  Section  IV. 

Fracture  Surface  Observation 

C-Ti  alloys  containing  57% 

Gamma  phase 

Figures  67  through  70  show  the  fracture  surface  of  C-Ti 
alloys  containing  57%  Gamma  phase  tested  at  different 
temperatures  ranging  from  25  degrees  C  to  550  degrees  C.  The 
arrow  shown  on  each  photograph  indicates  the  slip  direction 
during  the  compression  test.  Examination  of  the  fractographs 
indicate  that  the  amount  of  cleavage  failure  decreases  with  the 
increasing  temperature.  There  is  a  significant  amount  of  smooth 
slip  surface  in  specimens  which  were  tested  at  temperature  above 
450  degrees  C.  There  is  an  abundance  of  secondary  cracks  which 
indicates  that  significant  energy  was  expended  in  the  initiation 
and  propagation  of  these  cracks  prior  to  final  failure.  This 
phenomenon  existed  in  specimens  with  test  temperature  below  350 


degrees  C.  All  specimens  exhibited  mixed  mode  fracture 
(intergranular  and  transgranular) ;  but  the  dominant  fracture  mode 
was  intergranular.  These  features  indicate  that  failure  occured 
along  the  boundary  between  Alpha-Ti  and  Gamma  phase. 


C-Ti  alloys  containing  10% 

Gamma  phase 

The  fracture  features  of  C-Ti  alloys  containing  10%  Gamma 
phase  are  shown  in  Figures  71  through  74.  Only  a  small  amount  of 
cleavage  failure  and  secondary  cracking  are  observed  in  these 
specimens.  A  large  amount  of  smooth  slop  surface  exists  in  all 
specimens.  The  difference  between  specimens  below  and  above  350 
degrees  C  is  that  there  are  several  large  cracks  existing  in 
specimens  with  temperatures  below  350  degrees  C. 


Fractu-e  Features  of  C-Ti  Alloys 

Containing  10%  Gamma  Phase 

at  T-  350°  C  (A)  Complete  Cross  Section 

(B)  (C)  Higher  Magnification 


The  above  observations  occurred  at  elevated  temperature  as 
well  as  at  room  temperature.  Figure  60  and  61  display  that  the 
yield  stress  and  true  maximum  flow  stress  increases  with  volume 
fraction  of  Gamma  phase  at  constant  test  temperature.  Figure  62 
exihibits  that  the  true  failure  strain  drops  much  more  sharply  at 
temperatures  below  350  degrees  C  than  at  temperatures  above  350 
degrees  C.  The  fractographs  display  that  smooth  slip  surface 
exist  in  C-Ti  alloys  containing  low  volume  fraction  of  Gamma 
phase;  the  amount  of  cleavage  failure  largely  increases  in  C-Ti 
alloys  containing  high  volume  fraction  of  Gamma  phase  (test 
temperature  below  350  degrees  C) .  The  mobility  of  dislocation 
decreases  rapidly  as  test  temperature  decreases.  This 
information  indicates  that  ductility  of  C-Ti  alloys  with  high 
volume  fraction  of  Gamma  phase  decreases  abruptly  at  test 
temperature  below  350  degrees  C. 

Compression  test  results  display  that  the  true  maximum  flow 
stress  decreases  with  increasing  test  temperature  at  constant 
volume  fraction  of  Gamma  phase.  This  decrease  is  apprently  due 
to  the  annealing  and  recovery  of  the  dislocation  structure  at 
elevated  temperature  resulting  in  a  decreased  strain  hardening 
ability  of  the  material.  A  thermally  excited  lattice  structure 
assists  deformation  and  reduces  strength  at  elevated 
temperature.  The  fractographs  show  that  the  amount  of  secondary 
cracks  increases  as  test  temperature  decreases.  The  amount  of 
secondary  cracks  that  exist  is  indicative  of  the  amount  of  energy 


expended  in  the  initiation  and  propagation  of  these  cracks  prior 
to  failure.  Thus,  the  yield  stress  increases  as  test  temperature 
decreases  at  constant  volume  fraction  of  Gamma  phase.  The  amount 
of  smooth  slip  surface  increases  with  increasing  test 
temperature.  This  information  indicates  that  the  mobility  of 
dislocations  increases  as  test  temperature  increases.  This 
phenomenon  causes  the  true  failure  strain  to  decrease  with 
decreasing  test  temperature. 

Figure  64  (Y.S.  vs  1/T)  denotes  a  dependence  of  the  yield 
stress  (oy)  on  temperature.  A  phenomenological  relation  commonly 
used  to  establish  a  law  governing  the  deformation  behavior  is  the 
dependence  of  the  yield  stress  on  temperature,  which  can  be 
expressed  as 

oy  -  A  exp(Q/RT)  . (73) 

Where  A  is  a  constant,  R  the  Universal  gas  constant,  Q  an 
activation  energy,  and  T  the  absolute  temperature.  A  similar 
observation  was  applied  to  polycrystalline  TiC  by  Lipsitt  [68], 
Changes  in  slope  were  observed  in  plots  of  true  maximum  flow 
stress,  yield  stress,  and  true  failure  strain  vs  1/T. 
Re-examination  of  previous  data  [60,68,70]  from  both 
single-crystal  and  polycrystalline  TiC  hint  the  existence  of  such 
a  change  in  slope.  The  observed  changes  in  slope  suggest  that 
there  is  a  change  in  the  rate-controlling  mechanism  of 
deformation  at  and  above  the  transition  temperature  (T.T.).  T.T. 
was  found  to  increase  linearly  with  increasing  volume  fraction  of 


Gamma  phase.  An  increase  in  the  volume  fraction  of  Gamma  phase 
increases  the  number  of  carbon  atoms  that  will  become 
interstitial  atoms  in  a  C-Ti  alloy.  Ogden  [58]  indicates  that 
the  carbon  interstitials  in  solid  solution  in  titanium  cause 
titanium  to  undergo  a  brittle-to-ductitle  transition.  The  higher 
the  carbon  interstitials  content,  the  higher  is  the  transition 
temperature.  The  mechanism  which  causes  the  transition 
temperature  to  occur  in  our  study  is  still  unclear.  It  is 
suggested  that  transmission  electron  microscopy  be  used  to 
examine  regions  of  deformation  in  the  Gamma  particles  and  the 
grain  boundaries  to  help  define  the  deformation  mechanisms  and 
improve  the  properties  of  C-Ti  alloys.  The  results  of  this  study 
suggest  the  mechanical  properties  of  C-TI  alloys  are  controlled 
by  two  factors,  the  volume  fraction  of  Gamma  phase  and  test 
temperature . 

To  examine  the  effect  of  Gamma  phase  size  on  the  observed 
properties,  the  Hall-Petch  relation  (<7y  -  a0  +  d'1/^)  was  used. 

A  similar  observation  was  reported  by  Lipsitt  [68]  that  the 
dependence  of  the  yield  stress  on  grain  size  followed  a 
Hall-Petch  type  relation.  This  relationship  was  described  in 
Section  III.  A  plot  of  <7y/ d*1/2  vs  Gamma  phase  fraction  is  shown 
in  Figure  75.  By  using  a  statistical  program,  the  mathematic 
relation  between  yield  stress,  grain  size,  and  volume  fraction  of 
Gamma  phase  was  developed.  This  relation  can  be  expressed  as 


ay  -  A(d-1/2)*[exp(B*XP)] 


(74) 
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Figure  75.  Y.S./d*1/^  as  a  Function  of  Volume  Fraction  of  Gamma  Phase 


where  A  and  B  are  constants,  Xr  the  volume  fraction  of  Gamma 
phase.  Figure  76  is  the  plot  of  Cy/d'1/2  as  a  function  of  volume 
fraction  of  Gamma  phase  and  test  temperature. 

A  unified  equation  for  the  dependence  of  the  yield  stress  on 
grain  size,  Gamma  phase  fraction,  and  test  temperature  was 
developed  from  equations  (73)  and  (74).  The  relation  can  be 
expressed  as 

oy  -  A  (d-1/2)*[exp(B*Xp)]*[exp(Q/RT)] .  ---(75) 

Equation  (75)  indicates  that  we  can  obtain  the  yield  stress  of 
C-Ti  alloys  by  controllong  grain  size,  volume  fraction  of  Gamma 
phase,  and  test  temperature. 

Heat  treatment  studies  of  the  alloys  is  suggested  to 
determine  whether  inhomogeneities  can  be  eliminated  and  a  single 
phase  Gamma  alloy  produced  with  improved  ductility.  By 
controlling  grain  size,  volume  fraction  of  Gamma  phase,  test 
temperature  and  reducing  the  inhomogeneities  and  porosity  it  is 
envisioned  that  a  totally  single  phase,  high  carbon  content  off 
stoichiometric  material  can  be  produced  with  desired  high 
temperature  properties  and  ductility. 


Function  of  %Gamma 


I IV.  CONCLUSION 


The  following  conclusions  were  drawn  from  the  experiment  and 
analysis : 

0  The  yield  stress  increases  with  increasing  volume  fraction 
of  Gamma  phase . 

0  The  true  failure  strain  decreases  as  volume  fraction  of 
Gamma  phase  increases. 

0  The  true  maximum  flow  stress  and  yield  stress  decreases  as 
the  test  temperature  increases  at  constant  volume  fraction 
of  Gamma  phase . 

0  The  true  failure  strain  decreases  with  decreasing  test 
temperature  at  constant  volume  fraction  of  Gamma  phase. 

0  A  transition  in  the  controlling  deformation  mechanism  was 
found  to  occur  over  the  temperatuer  range  of  460° C  to  525° C. 

0  The  transition  temperature  increases  linearly  with 
increasing  volume  fraction  of  Gamma  phase. 

0  Fracture  surfaces  shows  a  mixed  mode  of  failure.  Regions  of 
cleavage  appear  to  be  the  result  of  failure  at  the 
Gamma-Alpha  phase  boundary. 

0  A  relation  between  yield  stress  (ay),  grain  size  (d) ,  volume 
fraction  of  Gamma  phase  (Xr) ,  and  test  temperature  (T)  was 
developed  to  describe  the  behavior  of  all  alloys  tested 
cr¥  -  A  (d'1/2)*[exp(B*Xr)  ]*[exp(Q/RT)  ] 
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THE  PARTICLE  SIZE  OF  GRMMR  PHRSE 
SPECIMEN  :  II 
TEMPERATURE:  25  deg.  C 
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THE  PARTICLE  SIZE  OF  GRMMR  PHRSE  DISTRIBUTION 
SPECIMEN  :  16 
TEMPERATURE:  25  deg.  C 
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THE  PHRTICLE  SIZE 
SPECIMEN  s  18 
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THE  PRRTICLE  SIZE  OF  GRMMR  PHRSE  DISTRIBUTION 
SPECIMEN  :  110 
TEMPERATURE:  25  dag.  C 


THE  PARTICLE  SIZE  OF  GRMMR  PHRSE  DISTRIBUTION 
SPECIMEN  s  112 
TEMPERATURE:  25  deg.  C 
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THE  PRRTICLE  SIZE  OF  GRMMR  PHRSE  DISTRIBUTION 
SPECIMEN  s  119 
TEMPERATURE:  25  deg.  C 


154 


evi 

< 

6 

£ 


X 

U 

x 

<r 


z 

o 

M 

)- 

D 

03 

*-» 

X 

h" 

tn 


u 

w 

x 

X 

D. 

X 

X 

r 

x 

o 

li. 

o 


o 


OJ 

o 

•a 


U  <S 
N  ru  in 
*-*  *-*  ru 

tn 


_J  u 

u  x 

w  z) 

h-  Z  H 
X  U  X 

x  r  x 

“■sis; 

ui  u  r 
X  X  Ul 
i-  o!  i- 


Aiznzausodd 


(SvUJUl)  866S2I 900  *  :b3dti  JO  30Ud3Ab 
BSE  s  *0N  313ZiaUd  IblOl 

*  *E’V9  *  NOZlDUdJ  3SbHd  U-«Mdlb 

*  99‘SE  J  NOZiDbdJ  3"!3Zidbd 


>  -V  u>. 


V' 

£ 


'  *1 


r .  ■ 

r: 


* 


u 


6 


*.1 


155 


3 


t  ^ 

I8- 

i  . 

V  3 


¥  3 


f  = 

I  I  ‘ 


5 


i 


-L_il_i- 


Ainiayaoyd 


(  g-vUim )  E8<iE98V00  *  :y3dy  JO  SDdiOAd 


662 

9i‘E9 

*2*9E 


*  on  3ioiiyyd  ibioi 

NOIlDtidJ  3SHHd  U-«Mdld 

Nonoyyj  3ioiidyd 


^.V-V-TwlV.. 


v,' 


-■-*  -  -  ■  - 


2 

o 


I- 

X 

m 


x 

h” 

in 


LJ 

in 

x 

x 

a. 


nj 

< 

E 

E 


X 

X 

X 

X 

u 


u 


X 

UJ 

X 

X 


o 


UJ 

N 


a> 

o 

TJ 


t-t  CM  If) 

tn  n  r\j 


u 

_J 

u 


UJ 

X 

X 


X 

X 

UJ 

X 

X 

UJ 


3 


t 


14 


;• 

i 


.1 


1 


.  «v 


i£v 


«•  (V  £ 

-  i  £ 

■  E  £ 

*  e  a: 

w  o  o 


<L  U. 

Ui  o  OJ 

O'  o 

a:  ui  -a 

n  <u 

M  04  in 

0)hn 


H  Z  H 
a:  u  <r 
trzo; 

OL  *-•  UJ 

u  a. 

UlilE 

IOLU 

h  cn  i~ 


Aimsusoad 


(g^uiui)  2(iSE8I900-  :B3BB  JO  33BB3AB 


89fr 

y.  ES  *  8fr 

y.  ^-ts 


’ON  33DIiyBd  IBiOl 
N0I1DBBJ  3SBHd  H-^MdlB 
N0I10BBJ  31DIlBBd 


k 


i1 


H 


THE  PRRTICLE  SIZE  OF  GRHMR  PHASE  DISTRIBUTION 
SPECIMEN  :  125 
TEMPERATURE:  25  dog.  C 
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SPECIMEN  :  126 
TEMPERATURE:  25  deg.  C 
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THE  PARTICLE  SIZE  OF  GAMMA  PHASE 
SPECIMEN  :  HA2 
TEMPERATURE:  150  deg.  C 
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TEMPERATURE:  150  deq.  C 
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THE  PARTICLE  SIZE  OF  GAMMA  PHASE  DISTRIBUTION 
SPECIMEN  :  HB2 
TEMPERATURE:  350  deg.  C 
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TEMPERATURE:  450  deg. 


PERRTURE:  450  deg. 
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THL  PRRTICLt  SIZE  OF  GRMMR 
SPECIMEN  :  HC6 
TEMPERRUJRE:  450  deq.  C 


SPECIMEN  :  HC? 
TEMPERATURE :  450  deg. 


TEMPERATURE:  550  deg. 
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THE  PARTICLE  SIZE  OF  GAMMA  PHASE 
SPECIMEN  :  H03 
TEMPERATURE:  550  deg.  C 


THE  PARTICLE  SIZE  OF  GRMMA  PHASE  DISTRIBUTION 
SPECIMEN  :  HD4 
TEMPERATURE:  550  deg.  C 
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SPECIMEN  :  HD5 
TEMPERATURE :  550  deg. 
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THE  PARTICLE  SIZE  OF  GRMMA  PHASE  DISTRIBUTION 
SPECIMEN  :  HDG 
TEMPERATURE:  550  deg.  C 


SPECIMEN  :  HD7 
TEMPERATURE:  550  deg. 


THE  PARTICLE  SIZE  OF  GRHMR  PHRSE  DISTRIBUTION 
SPECIMEN  :  HE 2 
TEMPERATURE:  B50  dog.  C 
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THE  PARTICLE  SIZE  OF  GRMMR  PHRSE  DISTRIBUTION 
SPECIMEN  :  HE5 
TEMPERATURE:  850  d«g.  C 
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